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We find an interpretation of the recent finding that the central charge density of the 
neutron is negative by using models of generalized parton distributions at zero skewness 
to relate the behavior of deep inelastic scattering quark distributions, evaluated at high 
X, to the transverse charge density evaluated at small distances. The key physical input 
of these models is the Drell-Yan-West relation We find that the d quarks dominate the 
neutron structure function for large values of Bjorken x, where the large longitudinal 
momentum of the struck quark has a significant impact on determining the center-of- 
momentum of the system, and thus the "center" of the nucleon in the transverse position 
plane. 

1. Outline 

Electron scattering is the preferred tool for extracting information on the spatial 
and momentum distribution of the quarks in nucleons. High energy scattering pro- 
vides a clean picture of the quarks' momentum distribution for a nucleon boosted 
into infinite momentum frame (IMF). Measurements of elastic scattering at lower 
energy scales allow extraction of the nucleon form factors, which can be related to 
the spatial distribution of charge in the rest frame of the nucleon. However, ob- 
taining rest frame charge distributions requires model-dependent relativistic boost 
corrections, limiting our ability to extract these distributions from data. Significant 
work has gone into better understanding the issues involved in studying nucleon 
distributions, as well as providing unified descriptions of the space and momentum 
distributions of the quarks. 

The present discussion and the papers on which it is based would not have been 
possible without the great amount of experimental technique, effort and ingenuity 
that has been used recently to measure the electromagnetic form factors of the 
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2 MILLER and ARRINGTON 

These quantities are probabihty amphtudes that the nucleon can 
absorb a given amount of momentum and remain in the ground state, and are 
related to the nucleon charge and magnetization densities. 

We note that there will be an Institute for Nuclear Theory Program 
held in the Fall of 2009 that is devoted to the electromagnetic physics 
of the Jefferson Laboratory upgrade to 12 GeV. Please see the website: 
www.int.washington.edu/PROGRAMS/09-03.html. 

We begin our analysis by reviewing recent worlj^l that determines the transverse 
charge density of the neutron in a model independent way. This work surprisingly 
found that the central charge density is negative. Then we discuss how the inclusive- 
exclusive connection is used to provide an interpretation for this facP. 

2. Definitions and the density interpretation 

We begin by by presenting definitions of the form factors. Let J^^{x'^) be the elec- 
tromagnetic current operator, in units of the proton charge. Then the nucleon form 
factors are given by the matrix element 

(p',A'|J''(0)b,A) =S(p',A') (^^^^F,{Q^)^i^qMQ^))u{p,\), (1) 

where M is the nucleon mass, and the momentum transfer qa = p'a— Pa is taken 
as space-like, so that = —q^ > 0. The nucleon polarization states are chosen 
to be those of definite light-cone helicities A, A'.'^ The charge (Dirac) form factor 
is Fi, normalized such that Fi(0) is the nucleon charge, and the magnetic (Pauli) 
form factor is F2, normalized such that -F'2(0) is the anomalous magnetic moment. 
The Sachs form factors are GsiQ^) = Fi{Q^) - ^F2{Q^), Gm{Q^) = Fi{Q^) + 

In the Breit frame, in which p = — p', Ge is the nucleon helicity flip matrix ele- 
ment of J°. This has been interpreted as meaning that G^; is the three-dimensional 
Fourier transform of the charge density in the rest frame. Indeed, the scattering of 
neutrons from the electron cloud of atoms measures the derivative —^dGE{Q^)/dQ^ 
at — 0, which has been widely interpreted as the mean-square charge radius of 
the neutron. However, a direct probability or density interpretation of Gb is spoiled 
by a non-zero value of Q^, no matter how small. This is because relativistic dy- 
namics, which cause the wave functions of the initial and final nuclcons of different 
momenta to differ, must be used. The final wave function is related to the initial 
one by a complicated boost operator that when acting on the initial state of a given 
momentum changes it to the same of the different final momentum. In general the 
boost operator (or boost) contains the effects of interactions. Thus the initial and 
final states differ, invalidating a probability or density interpretation. It is only in 
the case that non-relativistic dynamics are applicable that form factors are simply 
the Fourier transforms of the rest frame spatial distributions. 

It has commonly been assumed that at low momentum transfers, these cor- 
rections could be safely neglected. If we treat relativistic corrections as in atomic 
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physics, they are governed by v'^/c^ ~ /m,^, where m is the mass of the boosted 
constituent. Computing the form factors involves integrating over the momentum 
p so that it is replaced by the momentum transfer Q. Thus for small values of Q"^ 
we find 



where C is an unknown coefficient. In a simple constituent quark model, m « 
300 McV, and so these 1/m^ boost corrections are negligible only forp^ <C 0.1 GcV^. 
However, for the neutron Ge{0) = 0, so the corrections are expected to be rela- 
tively large. The boost correction term need not be small compared to the finite 
size contribution unless C <C 1; since both corrections scale with Q^, the relative 
correction to the extracted radius does not vanish as approaches 0. The con- 
stituent quarks represent the low Q^, dressed versions of the near-massless current 
quarks of QCD. For current quarks of mass 5-10 MeV, the boost corrections are 
important for all values where measurements exist. While it is possible to con- 
struct systematically improveable models which should yield a complete description 
of the nucleon, it is not clear how one would quantitatively determine how well any 
particular model fully reproduces the corrections associated with the rolativistic 
boost. Thus, there is always some model dependence in the extraction of the rest 
frame charge distributions from the form factors, and it is not clear how well one 
can quantiiy these corrections and uncertainties, even at very low Q"^. 

3. Light Cone Coordinates, the Infinite Momentum Frame, and 
the Drell-Yan Frame 

The use of light cone coordinates and the kinematic subgroup of the Poincare group, 
which is closely related to the use of the infinite momentum frame, enables one to 
avoid the difficulties associated with including the boost. 

The basic idea is that the "time" variable is given by a;"*" = {ct + z)l\f2 = 
[x^ + x'^)/ a/2 and the conjugate evolution operator is given by p~ = (p° — p^) / \/2. 
Starting in an ordinary reference frame and making a Lorentz transformation into 
a frame moving with nearly the speed of light in the 3 direction converts the usual 
t into x'^ . The 3 spatial variables must be different than so we take x~ = 
{x^ — x^)/\/2. If a;"*" = then x~ = —\f2z and x~ can be thought of as something 
like a z or a;^ variable, but rotational invariance can not generally be used to relate 
the and .t, y dependence of the density. The canonically conjugate momentum 
to x" is p+ = -|-p^)/a/2. For the transverse degrees of freedom we use the usual 
position (cc, y) and momentum {px,Py) variables denoted by b, p where the boldface 
notation denotes a transverse vector. 

We also exploit the kinematic subgroup. That is transverse Lorentz transfor- 
mations, with a transverse velocity v (or boosts in the transverse x,y or b direc- 
tion) do not involve interactions. In particular, these transformations are defined 
by fc+ — * k — + k — k~^v, with k~ changed so that fc^fc'* is not changed. These 





March 9, 2009 18:20 WSPC/INSTRUCTION FILE pionijmpe 



4 MILLER and ARRINGTON 

transformations are just like the non-relativistic Gahlei transformation except that 
appears instead of a mass. This means that we are allowed to do Fourier trans- 
formations of variables provided transverse degrees of freedom are involved. 

If the momentum transfer vector is space-like we may use the so-called Drell- 
Yan (DY) frame in which — 0. This means that the plus component of the 
momentum of the nucleon is the same before and after the absorption of the single 
photon. Then the momentum transfer vector is in the transverse direction, so that 
usual two-dimensional Fourier transform techniques can be used. 

4. Charge density of the neutron 

A proper determination of a charge density requires the measurement of a density 
operator. We shall show that measurements of the pion form factor directly involve 
the three-dimensional parton charge density operator, in the infinite momentum 
frame IMF, /5oo(a;^,b) = J+(a;^,b). In this frame the electromagnetic charge den- 
sity J° becomes J+ and 

Poo{x~,h) ^^eqq(x^ ,h)-i+ q{x^ ,h) = ^ eqV2q^(x", b)(7+(x", b), (3) 
q 9 

where q+{x^) = 7°7^/\/2<z(a;'^), the independent part of the quark-field operator 
q{x^). We set the time variable, x'^ = [t + z)/\/2, to zero, and do not display it in 
any function. ,„„^ 

The spatial structure of a hadron can be examined if one useJ^^^H^ states that 
are transversely localized. The state with transverse center of mass R set to is 
formed by taking a linear superposition of states of transverse momentum: 

\n = Q,\)^M j ^|p+,p,A), (4) 

where , p. A) are light-cone helicity eigenstates and A/" is a normalization factor 
satisfying \}J'\^ J ^ ^' expansion Eq. ^ makes sense only in the infinite 
momentum frame — > oo because we must have 2p^p^ = M'^. The nucleon 
states are normalized as (y*", p'A'|p+, p, A) = 2p+ {2'k)^5x' ^x6{p'^ - p'^)S'^^^p' - p) . 
Next we relate the charge density 

R = 0, A| p^{x-,h) b+, R = 0, A) 
'^'^ (p+,R = 0,Ab+,R = 0,A) ' 

to Fi(Q^). In the DY frame no momentum is transferred in the plus-direction, so 
that information regarding the x~ dependence of the distribution is not accessible. 
Therefore we integrate over x~ , using the relationship translational invariance, and 
then use our momentum expansion and the definition of Fi as the helicity non-flip 
matrix element of J+ in a Drell-Yan frame to find 



\P 



p{b) ^ J dx poo{x ,6) (6) 
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Fig. 1. Upper panel: proton transverse charge density p{b). Lower panel: neutron transverse charge 
density. The solid curves use the parameterization of Kelly, and the dashed (red) curves use 
Bradford, et ah. Reprinted with permission from Ref. ^.Copyright 2007 by the American Physical 
Society. 

where Jq is a cylindrical Bessel function and p{b) is termed the transverse charge 
density, giving the charge density at a transverse position b, integrating over aU lon- 
gitudinal momentum. In the infinite momentum frame, the longitudinal dimension 
of the nuclcon is contracted to a point, and only the transverse position remains. 
The value of 6 = corresponds to the center of longitudinal-momentum of the 
nucleon in the transverse dimension. 

It is not as straightforward to isolate t he rnag netization distribution, and there 
have been multiple such densities proposed^UGST Starting with the matrix element 
of J • A and working in the infinite momentum frame leads to the result that the 
transverse magnetization density to is the two-dim ens ional Fourier transform of ^2, 
just as the charge density is the transform of Fi This interpretation yields a 
difference between the magnetic and electric radii in the prot oii. 

Here we exploit Eq. ([7]) by using recent parameterizations^^^^I^Iiniof measured 
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form factors to determine p(b ). App lying Eq. ([7]) to the proton using two sets of 
form factor parameterizationJl^IEl yields the results shown in the upper panel 
of Fig. [1] The curves obtained using the two different parameterizations overlap. 
Furthermore, there is negligible sensitivity to form factors at very high values of 
that are currently unmeasured. The density is peaked at low values of 6, but has a 
long positive tail, suggestive of a long-ranged, positively charged pion cloud. 

The neutron results for p{b) are shown in the lower panel of Fig. [1] The curves 
obtained using the two different parameterizations seem to overlap, but we will 
return to this below. The surprising result is that the central neutron charge density 
is negative. The negative nature of the neutron's central charge density appears to 
contradict two current ideas. If the neutron is sometimes a proton surrounded by 
a nega tive ly charged pionic cloud, one would expect to obtain a positive central 
densitjff^. Another mechanism involving correlations in the nucleonic wave function 
induced by one gluon exchange would also lead to a positive central density because 
the interaction between two identical d quarkJI^ is repulsive. 



The resultant negative central density thus deserves further examination. The 
upper panel of F ig. s hows Fi for the neutron obtained using the two different 
parameterizationsi4113|-yirhich are observably different. However, in both cases. -Fi is 
negative for all values of Q^. If Fi is always negative, then taking 6 = 0, Jo{Qb) = 1 
in Eq. ([7]) , will always yield a negative central density. The long range structure of 
the charge density is captured by displaying the quantity bp(b) in the lower panel 
of Fig. [2] At very large distances from the center, the charge density is negative, as 
expected in the pion cloud picture. 

These findings appear to contradict previous understanding of the nucleon 
charge distributions based on the model-dependent extraction of the rest frame 
charge distributions. Because of the large and model-dependent boost corrections 
at large Q^, one cannot obtain information about the rest frame charge density at 
the very center of the nucleon from measurements of the form factors. While there is 
no direct experimental information on this possibility of a very small negative core 
in the rest frame distribution, this negative core seems to contradict the accepted 
explanations of the origin of the long range negative cloud which occurs along with 
a positively charged interior. In addition, the negative core in the IMF transverse 
density is a feature even in models that build up the charge distribution based on 
a pion cloud model. It is clearly important to understand the differences between 
the IMF charge density and the rest frame charge density to fully understand the 
new features of these model-independent spatial distributions. 

The surprising model independent result is that the density of the neutron is 
negative. The remainder of this presentation is concerned with trying to explain 
this remarkable feature of nature. 
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Fig. 2. Neutron Fi and transverse charge density. Upper panel: _Fi. Lower panel: bp{b) in transverse 
position space. The solid curves are obtained using the fits of Kelly, and the dashed curves the 
fits of Bradford, et al.. Reprinted with permission from Ref. ®. Copyright 2007 by the American 
Physical Society. 



5. Inclusive-exclusive connection 

Generalized parton distributions (GPDs) contain information about the longitudi- 
nal momentum fraction x as well as the transverse position b. Information regarding 
the X and h dependence is obtained from experiment by using GPDs to reproduce 
both deep inelastic scattering and elastic scattering data. Thus we use this inclusive- 
exclusive connection to bette r unde rstand the central neutron charge density. 

The widely studied are of high current interest because they can 

be related to the total angular momentum carried by quarks in the nucleon. We 
consider the specific case in which the longitudinal momentum transfer ^ is zero, 
and the initial and final nucleon helicities are identical (A' = A) . Then, in the light- 
cone gauge, = 0, the matrix element defining the GPD Hq for a quark of flavor 
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q and zero skewness is 

Hgix,t) = {p+,p',X\dgix,0)\p+,p,X), (8) 

where 

0,(x,b) ^ J (-^,b) g+ (^,b) e"^"--. (9) 

We abbreviate Hq{x, ^=0, t) = Hg{x, t) and -t = -(p' - pf = (p' - p)^ Q'^. 

It is w ell k nown that GPDs provide a unified description of a number of hadronic 
properties!!^ Of particular interest here is that for t=0 they reduce to conventional 
PDFs, Hq{x, 0) = q{x), and that the integration of the charge- weighted Hq over x 
yields the nucleon electromagnetic form factor: 

Flit) = y^eq [ dxHqix,t). (10) 



The impact parameter-dependent PDI^^QI for a quark of flavor q is the matrix 
element of the operator Oq in the state R = 0, A): 

plih, x) = R = 0, A| dq{x, b) R = 0, A) . (11) 

We use the notation (b, x) instead of the originally definecJ^Hl q^x, b) because this 
quantity is a density that gives the probability that the quark has a longitudinal 
momentum fraction x and is at a transverse position b. The quantity p''^{h,x) is 
the two-dimensional Fourier transform of the GPD Hq-. 

pl{h,x)= [ -fl-e~^^-^Hq{x,t^-<i'). (12) 



(2n) 



We extract the form factor Fi by integrating p^(b,a;) over all values of 



X, 



multiplying by the quark charge Cq, and summing over quark flavors g.'^'The result 
is the IMF charge density in transverse space: 

p(6)^^e, /'da.pKb,^). (13) 

The relations Eq. (|13p and Eq. ^ provide two expressions for the transverse den- 
sity p{b). The impact parameter GPD and the three-dimensional density are related 
by Parseval's theorem. The quantity p{b) gives the charge density at a transverse 
position b irrespective of the longitudinal momentum fraction or longitudinal po- 
sition. 

There is a tight connection between the values of x and the values of b. For a 
given Fock space component, the center of transverse momentum (CM) R is given 

by 

R = = ^x,b„ (14) 
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where the longitudinal momentum fractions and transverse positions of 

the i'th partons. Using this definition, the longitudinal momentum of a quark de- 
termines it's impact on defining the transverse position of the nucleon. If the struck 
quark, Xi , carries nearly all of the plus component of the total momentum, then the 
other quarks must carry a total plus momentum of 1 — Xi which approaches zero 
as — > 1. Thus, Xj — > for i ^ j, and via Eq. (|14|) . must also approach zero as 
long as the values of bj remain finite. Thus, large values of x correspond to a small 
value of 6, as the struck quark plays a large role in defining the transverse center 
of mass as X ^ 1. 

Let us investigate p_L(b,a::) to understand the origin of the neutron's negative 
central charge density. The quantities are not measured directly, but have been ob- 
tained from models t hat incorpor ate fits to parton distributions and electromagnetic 
nucleon form factors. I^IElESlMl'phis method exploits form factor sum rules at zero 
skewness to model information regarding the valence quark GPDs, = H"^ — H''. 
This yields the net contribution to the form factors from quarks and anti-quarks, 
although it does not correspond to the valence distribution within a model for which 
sea distributions for quarks and antiquarks have different x or i dependences. The 
effects of strangeness are neglect ed in thes e fits. 

Each parameterization we usJ^lESES] incorporates the Drell-Yan-WesiPSl rela- 
tionship between the behavior of the structure function i'W2 (x) function near x — I, 
measured in inclusive reactions and the behavior of the electromagnetic form fac- 
tor at large values of Q^, measured in the exclusive elastic scattering process. In 
particular, for a system of n + 1 valence quarks, described by a power-law wave 
function 

limz.T4/2(x) = (l-x)2"-i^ lim F^{Q') ^ (15) 

with the relation being that the same value of n that defines the high- a; behavior 
of the structure function also defines the high-Q'^ behavior of the form factor, thus 
associating the behavior of large values of x with large momentum transfers, — 
, which in turn correspond to small values of b. So again we see we see a connection 
between large values of x and small values of b. 

To proceed further we use specific forms of the GPDs, and these determine the 
details of the results. Diehl et aL^^H 

H^{x,t) = qy{x)exp[fq{x)t], (16) 

where 

f,{x) = [a'\og[l/x]+B,]il^xf +Ayx{l-x)^, (17) 

is the form that gives the best fit to the data. The parameter a' represents the slope 
of the Regge trajectory, and the CTEQ6 PDF*^ are taken as input. Here we use 
the best fit parameters, taken from the second line of Table 8 of Ref. ^221 These are 
detailed in Ref. The labels q refer to u and d, the u and d quarks in the proton. 
These correspond to d and u quarks in the neutron, if charge symmetry l ^7 | 28 | 29 |30j 
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is upheld. For the proton, Id^juy faUs rapidly for large values of x, which means that 
u quarks dominate the parton distribution for large values of x. For the neutron, the 
assumption of charge symmetry implies that for the neutron is the same as Uy 
for the proton, and vice- versa. Thus, the d quarks dominate t he p arton distribution 
in the neutron for large values of x. The distributions of Ref. have Aq = Bq = ^ 
and jqix) — [a^ log[l/a;]](l — x). Those of Ref. have a more complicated form 
and include the additional constraint that the nucleon consists of three quarks at 
an initial scale of Ql = 0.094 GeV^. 

We study the connection between regions of x and regions of h. To do this define 

pl{b,Ax)= [ dxeqpl{b,x), (18) 

J Ax 

where eq is the quark charge in units of the proton charge (e„ = 2/3, — —1/3) 
with p^" being obtained from appropriate sums of p^. This represents the con- 
tribution to the charge density from quarks in the x region defined by Ax, rather 
than the total density obtained by integrating over all x. 




b [fm] 

Fig. 3. The proton transverse charge density, p'^{b, Ax), for quarks in different Ax regions: 3;<0.15 
(sohd), 0.15<a;<0.3 (long-dash), 0.3<a;<0.5 (short-dash), and a;>0.5 (dotted). The curves have 
been normalized to unity at b = to emphasize the variation in width. Reprinted with permission 
from Ref. ^. Copyright 2008 by the American Physical Society. 



We have pointed out above that at large x, the struck quark plays a significant 
role in defining the transverse CM, so the distribution of high-x quarks becomes 
localized at small values of b. This is clearly visible in Fig. [31 which shows p^{b, Ax) 
for different bins in Ax. The curves have been scaled to yield unity at 6 = 0, to 
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emphasize the variation in width. The four Ax regions yield 58%, 25%, 14%, and 
3% of the total charge, with the largest contributions coining from the bins with 
the smallest values of x. For x ~ 0.1, the half-maximum width is 0.5 fm, while for 
X « 0.8, it is 0.12 fm. Thus the large x quarks (mainly u quarks in the proton) play 
an increasingly prominent role in the charge distribution at s mall values of b. The 
curves shown in Fig. [3] are obtained using the GPD of Ref. I^^l The results obtained 
from the Guidal et al, parameterization for the GPDs are barely distinguishable. 
The GPDs of Ref. ■ also have a strong tendency to be constrained to smaller and 
smaller values of b as the value of x increases. We evaluate the GPDs of all three 
models using the starting scale Qq of each model. 




Fig. 4. Ratio of u quarks to d quarks in the neutron from several analyses of deuteron and proton 
data. The solid line is the CTEQ6L parameterization 



Taking what we have learned from the proton, we now consider the neutron. 
Figure |4] shows the ratio of the up- to down-quark distributions in the ne utron, 
as extracted from various analyses of deuteron and proton data I 31 | 32 | 33 | (^gpg- 
different models for the nuclear corrections in dueterium), and from the CTEQ6lI^ 
parameterization. For a; ^ 0, the up and down quark distributions are similar, and 
the contribution to the charge distribution from this limit should be similar to 
that of the proton; a broad distribution of net positive charge. For x = 0.3 and 
above, the u quark distribution is less than half the d quark distribution, yielding 
a net negative contribution to the charge. Because the distribution of quarks is 
more localized near & = as x increases, a negative peak can be formed if there 
is a sufficiently large contribution from down quarks at large x values. Above x = 
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0.5, dn{x)/un{x) is at least three, except for the density-based extrapolation (blue 



circles), and continues to increase with x. As discussed in Ref. even if one 

uses a density-based extrapolation of the EMC effect to apply nuclear corrections, 

the implementation used in this analysis significantly overestimates the effect for 

deuterium, and yields an unrealistically large result at very large x values. In this 

region, the net impact on the charge distribution will be negative, and will be 

peaked at smaller values of b. This is shown in Fig. [5l which separately shows the 

contrib utio ns to the neutron charge de nsity fr om u and d quarks based on the GPD 

of Ref. xhe distributions of Refs. I ^-*^ ! ^*^ ! yield somewhat different results, but 

[991 

exhibit the same qualitative behavior. For example, the GPDs of Ref. shown 
in Fig. [5l yield a negative central neutron charge density for valu es o f x between 
0.15 and 0.3 and between 0.3 and 0.465, but for the GPDs of Ref. 1211, the central 
density is positive unless x is slightly greater than 0.465. 
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Fig. 5. The u and d quark contributions to the neutron transverse charge density, p" (6, Ax) and 
p^(6, Ax). Here the quark flavor refers to the neutron {u in the proton is d in the neutron). The 
curves correspond to the same Ax regions as in Fig. [3] The largest contributions come from small 
X, where u and d quarks contribute similar amounts of charge. As one goes to larger x values, the 
charge is shifted to smaller values of fe, while at the same time the up quark distribution drops 
rapidly with respect to the down quarks. Reprinted with permission from Ref. ^. Copyright 2008 
by the American Physical Society. 



Next we examine the total charge distribution of the neutron. Fig. [6] shows the 
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charge distribution of the neutron, separating out the contributions from low and 
high X regions, and shows fop" to suppress the very large density at the center. For 
X < 0.23 the charge distribution (dotted line) is positive for b < 1.5 fm and slightly 
negative distribution at larger radii. For x > 0.23, the contribution (dashed line) is 
largely negative, and highly localized below 0.5 fm. The negative region at the center 
of the neutron transverse charge distribution arises a natural consequence of the 
model-independent definition of the charge density. The low momentum partons 
have a larger spatial extent and reproduce the intuitive result of the pion cloud 
picture: a positive core with a small negative tail at large distances, although the 
negative tail is difficult to see for this parameterization of the GPD, given the large 
scale required to show the negative core at small b. 



T I I I I I I I I I I r 




J 1 I I I L 



0.0 0.5 1.0 1.5 

b [fm] 

Fig. 6. Transverse charge density for the neutron. The dotted Une is the contribution from x < 0.23, 
dashed is that for x > 0.23, and the sohd is the total. Reprinted with permission from Ref. ^. 
Copyright 2008 by the American Physical Society. 



We also study the quantity 

p^{b,x) = J2 ^.Plib,^) (19) 

to obtain a pictorial view of the transverse charge density for specific values of x 
(Fig. [7]). The striking feature of the negative spike appears prominently for x = 0.3 
and more prominently for x — 0.5. These figures show how the central negative 
charge density appears more and more prominent as x increases. Clearly the nega- 
tively charged d quarks dominate at the center of the neutron. 



March 9, 2009 18:20 WSPC/INSTRUCTION FILE pionijmpe 



14 MILLER and ARRINGTON 




Fig. 7. Neutron generalized parton distributions p{b,x) for x = 0.1 (top), x = 0.3 (middle), and 
X = 0.5 (bottom). 

As discussed earlier, the longitudinal-momentum weighting used in determining 
the nucleon center of mass leads to a strong correlation between the position of the 
struck quark and the center of the nucleon for large x. It is informative to try and 
remove this effect to obtain something that is closer to our intuitive picture. We 
can do this by examining the position of the struck quark relative to the center of 
the spectator system, so that the struck quark does not influence the definition of 
the center of the neutron. This can be approximated by looking at the position of 
the struck quark relative to the spectators. We use Eq. ([T4|) with the origin set to 
the center of momentum, for a struck quark at (xi,bi) = (a;,b), to determine the 
momentum- weighted spectator position, bgpcc, and the relative distance from the 
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struck quark to the spectator quarks: 

xibi + ^Xibi = xb + (1 - x)bspoc = 0, (20) 

i>l 

Brcl = b — bspcc = jz r = Brd- (21) 

(1 - x) 

We exhibit the dependence on Bid by defining a function 

p'^/^'iBrel^x) = p^{Brel{l - X),X) (22) 

which gives the probability that a struck quark of longitudinal momentum fraction x 
is a distance Brei away from the spectator center of momentum. Figure [5] shows this 
rescaled version of p±{b), with the contribution at each x value normalized to unity 
at & = 0. The quantity p^^°'^(Bi.ei, x) does not correspond to a true density, but can 
provide a better approximation to our intuitive picture of the charge distribution, 
as it removes the influence of the struck quark on defining the center of the nucleon. 
While the charge distribution coming from very low x quarks has a greater spatial 
extent, the decreasing width of the p± (b) distribution for large x quarks is essentially 
completely removed when looking at -Bid- 




Fig. 8. The u and d quark contributions to pj'""''"(B^ai,x) see Eq. I I22I I. vs B^-^i for x = 0.1 (solid), 
0.3 (long-dasii), 0.5 (short-dash), and 0.7 (dotted). The curves are scaled to unity at B^ei = 0. 
Here the quark flavor refers to the neutron {u in the proton is d in the neutron). Reprinted with 
permission from Ref. ®. Copyright 2008 by the American Physical Society. 
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We s ummariz e our findings with the statement that, using the model GPDs 
of Refs. | 22 | 21 | 23 |^ ^j^g dominance of the neutron's d quarks at high values of x 
leads to a negative contribution to the charge density which, due to the definition 
of b, becomes localized near the center of mass of the neutron. This localization 
does not appear when examined as a function of the position of the struck quark 
relative to the spectators, and is an consequence of the fact that quarks with a large 
longitudinal momentum play an important role in defining the transverse position 
of the neutron. 
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